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Abstract. Consider the Navier-Stokes flow past a rotating obstacle with a general time- 
dependent angular velocity and a time-dependent outflow condition at infinity - sometimes 
called an Oseen condition. By a suitable change of coordinates the problem is transformed 
to an non-autonomous problem with unbounded drift terms on a fixed exterior domain 
O C K d . It is shown that the solution to the linearized problem is governed by a strongly 
continuous evolution system {Tn(i, s)}t> s >o on for 1 < p < oo. Moreover, L p -L q 

smoothing properties and gradient estimates of Tn(t, s), < s < t, are obtained. These 
results are the key ingredients to show local in time existence of mild solutions to the full 
nonlinear problem for p > d and initial value in L p a (VL). 

1. Introduction 

In this paper we consider the flow of an incompressible, viscous fluid in the exterior of 
a rotating obstacle subject to an additional time-dependent outflow condition at infinity. 
Here the angular velocity of the obstacle and the outflow condition at infinity may depend 
on time and also the axis of rotation may change. The equations describing this problem 
are the Navier-Stokes equations in a time-dependent exterior domain with a prescribed 
velocity field at infinity. 

After rewriting the problem on a fixed exterior domain Q C M. d , we obtain an non- 
autonomous system of equations involving a family of time-dependent operators of the 
form 

A(t)u = P n (Au + (M(t)x + c{t)) - Vu - M{t)u) , t>0. (1.1) 
where Pr> denotes the Helmholtz projection from L p (Q) d onto the solenoidal space L^(Q) 
and M E C^flO, oo); R dxd ), c E C^flO, oo); R d ). The main difficulty in dealing with these 
operators arises since the term M(t)x- V has unbounded coefficients in the exterior domain 
Q. In particular, the lower order terms cannot be treated by classical perturbation theory 
for the Stokes operator. In the autonomous case M{t) = M and c(t) = such an operator 
was first considered by Hishida [14j and then later by Geissert, Heck, Hieber [8]. It is the 
aim of this paper to extend their result to the non-autonomous case. 

In the following let us briefly motivate our problem and let us show why it is interesting 
to study the non-autonomous case. For this purpose let O C M. d be a compact obstacle 
with smooth boundary and let fl := M, d \ O be the exterior of the obstacle. Furthermore, 
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let m G C 1 ([0, oo); M dxd ) be a matrix- valued function that describes the velocity of the 
obstacle. Then, the exterior of the rotated obstacle at time t > is represented by 
Q(t) := Q(t)Q where Q G C 1 ([0, oo),M. dxd ) solves the ordinary differential equation 

d t Q{t) = m{t)Q{t), t>0, 

Q(0) = Id. 

With a prescribed velocity field G C 1 ([0, oo); M d ) at infinity, the equations for the fluid 
on the time-dependent domain Q(t) with no-slip boundary condition take the form 

v t - Av + v-Vv + Vq = in (0,oo) x Q(t), 

divv = in (0, oo) x f2(t), 

v(t, y) = m{t)y on (0, oo) x dfi(t), (1.3) 

lim v(t, y) = foo(^) for t G (0, oo), 

|y|-k>o 

v(0,y) = v (y) in SI. 

Here v = v(y,t) and q = q(y,t) are the unknown velocity field and the pressure of the 
fluid, respectively. 

As usual, it is reasonable to reduce ( II. 3p to a new problem on a fixed exterior domain by 
some suitable coordinate transformation. Since m(-) is the velocity of the rotated obstacle, 
it is natural to assume that m(t) is skew symmetric for all t > 0. This implies that for all 
t > the matrix Q(t) is orthogonal. Thus, we can set 

x = Q(t) T y, u(t,x) = Q(t) T (v(t,y)-v oc (t)), p(t, x) = q(t, y). (1.4) 

Then we obtain the following new system of equations on the reference domain Q: 

u t -Au- M (t)x ■ Vu + M(t)u 1 

+Q(t) T v OQ (t) ■ Vn + Q(t) T ^oo(t) 1=0 in(0,oo)xfi, 
+u ■ V« + Vp J 

div u = in (0, oo) x Q, 

u(t, x) = M{t)x - Qiffv^it) on (0, oo) x dQ, (1.5) 

lim u(t,x) = for t G (0,oo), 

\x\— >co 

u(0, x) = u (x) in Q. 

Here M(t) := Q(t) T m(t)Q(t) is the transformed velocity of the obstacle. The coordinate 
transformation also ensures that the new velocity field u vanishes at infinity, which is a 
natural condition in the L p -setting. 

Note that a problem of this type also arises in the analysis of a rotating body with a 
translational velocity — v^t) by a similar coordinate transformation, see e.g. the explana- 
tions in [3]. 

Problem (ll.3p was studied intensively for the special case of time-independent matrices 
M(t) = M and without an outflow condition, i.e. = 0. Hishida [H] showed that the 
solution to the linearized problem is governed by a strongly continuous semigroup on L 2 , 
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which is however not analytic. Moreover, he constructed local mild solutions in L 2 by using 
the Fujita-Kato approach (cf. jl]). Later this generation and existence result was extended 
to the general L p -theory by Geissert, Heck, Hieber [8]. Hishida and Shibata [16] were even 
able to show global extistence for small data. The model problem in M d was studied by 
Hieber and Sawada [13j in the L p -setting. 

The case of time-dependent angular velocities was considered by Hishida |T5] in the L 2 - 
context, however he assumes that the axis of rotation is fixed and that the sign of his 
angular velocity does not change. 

For the problem including an additional outflow condition at infinity, there are only a 
few results. The case, where M(t)x = u)(t) x x and uj : [0, oo) — > M 3 is the angular velocity 
of the obstacle and : [0, oo) — > M 3 a time-dependent outflow velocity was considered by 
Borchers [2] in the framework of weak solutions. This work was somehow the starting point 
in the analysis of viscous fluid flow past a rotating obstacle. More recently, Shibata [20] 
studied the special case where M(t) = M, Woo(t) = and Mv^ = 0. The additional 
condition Mv^ = 0, i.e. Q(t) T t>oo = kvoo f° r k G {+1,-1}, ensures that (jl.5p is still an 
autonomous equation. The physical meaning of the additional condition is that the outflow 
direction of the fluid is parallel to the axis of rotation of the obstacle. The stationary 
problem of this latter situation was analysed by Farwig [3] for the whole space case M 3 . 

In order to relax the assumption Mvoo = 0, i.e. in order to allow a general Voo, it is 
necessary to study a non-autonomous problem. Thus, in this context it is natural even to 
allow time-dependent outflow velocities i>oo(0 an d time-dependent angular velocities M(-). 
The non-autonomous model problem of (jl.5p in the case Q = M. d was recently studied by 
the first author [TT| and by Geissert and the first author [B]. Indeed, they were able to 
show that the family of operators in (II. ip . equipped with suitable domains, generate a 
strongly continuous evolution system on L^(M. d ), 1 < p < oo, which enjoys nice regularity 
properties. Their approach is based on an explicit solution formula for the linearized 
problem. By a version of Kato's iteration scheme (cf. OUT]) one obtains a (local) mild 
solution to the nonlinear problem on M. d for initial value Vq G L p r (M d ), d < p < oo. In this 
paper we use their results for the linearized problem to cover the physically more realistic 
situation of exterior domains by some cut-off techniques. 

Notations. The euclidian norm of x G M. d will be denoted by \x\. By B(R) we denote the 
open ball in M. d with centre at the origin and radius R. For T > we use the notations: 



Let us come to notation for function spaces. For a C ' domain O C 1 and 1 < p < 
oo, j G N, W^ p {Vi) denotes the classical Sobolev space of all L p (f2)-functions having weak 
derivatives in L P (Q) up to the order j. Its usual norm is denoted by || • \\j, p ,n and by || • \\ P) n 
when j = 0. If Q = M. d we drop Q in the notations of the above norms. We will use also 
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the notation 

(/,<?)q:= / fgdx, feL p (Q), g eL p '(n) 

with i + i = 1. By Wq' p (Q), s > 0, we denote the closure of the space of test functions 
C£°(f2) with respect to the norm of W S,P (Q). For s < we set 

w s ' p (fi) : = (w- s ' p \n) 

where - + = 1. We denote by H S ' P (Q) with s G (0, 2) the Bessel potential spaces, which 
are defined by complex interpolation 

H S ' P (Q) := [L p (Q),W 2 ' p (Q)] i . 

Its norm will be denoted by || • \\ s ,p,n- Moreover, we set 

C£(fl) := {/ G C?(n) d : div/ = 0}, 

L p (fi) d :={/ = (/i,..., / d ) :/ < GL ? (fi) ) i = l 1 ..., d}, := C~ (fi) , 

:= {Vp : p G W 1 '^)}, W^fi) := {p G Lf oc (H) : Vp G 

It is well-known that for a Lipschitz domain Ocl* 1 with compact boundary the Helmholtz 
decomposition holds (see e.g. j5] for more information): 

L p (Q) d = L p a (tt) © G p (n). 

The projection from L p (Q) d onto L p (Q) is denoted by Pq. 

2. Main results and strategy of proofs 

In this section we present the main results and sketch the basic strategy of the proofs. 
In the following O C M. d is always a compact obstacle with C 1,1 -boundary and Q := M. d \ O 
is an exterior domain. Moreover, M G C 1 ([0, oo), IR dxd ) and G C 1 ([0, oo), ~R d ) are as 
described in the Introduction. In particular recall that tr M(t) = for all t > 0. To 
simplify our notation, we set c(t) := -Qitfvooit). Since the term Q(t)dtVoo(t) in equation 
(ll.Sp is constant in space, we may put this term in the pressure p. Thus, in the following 
we consider the system 

u t - Au - (M(t)x + c{t)) ■ Vw + M{t)u 
+u-Vu + Vp 



= in (0, oo) x Q, 

div u = in (0, oo) x Q, 

u(t,x) = M(t)x + c(t) on (0, oo) x <9fi, (2.1) 

lim u(t,x) = for t G (0, oo), 

\x\— ¥oo 

u(Q, x) = uo(x) in Q, 

where divu = 0. Moreover, we assume that the initial value Uq satisfies the compatibility 

assumption Uq ■ v = (M(0)x + c(0)) • v on dfl, where v denotes the outer normal vector. 
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As a first step we construct a solenoidal extension in Q of the boundary velocity u(t, x)\qq. 
For this purpose we introduce the Bogovskii operator, which concerns the solution of the 
equation div u = f in appropriate function spaces. This operator will also be needed later 
in Section [5] to keep the solenoidal condition in our cut-off procedure. For proofs and more 
information on the Bogovskii operator we refer to [TJ[7] and to the monograph [5]. 

Lemma 2.1. Let D C M> d , d > 2, be a bounded Lipschitz domain, 1 < p < oo and k G No- 

(a) There exists a continuous operator 

M D : Wq' p (D) (Wo +1 ' p (D)) d 

such that 

divB D / = / 

for all f G Wq' p (D) satisfying f D fdx = 0. 

(b) For k > —2 + -, the above operator IB o can be continuously extended to a bounded 

operator from W^ P (D) to (W" fe+1 ' ? \D)) d . 

Let C G C^°(M. d ) be a cut-off function with < ( < 1 and £ = 1 near dQ. Moreover we 
set K := supp VC and define b : [0, oo) x R d — y R d by 

b{t, x) := ({x){M{t)x + c(t)) - B^((VC) ■ {M{t)x + c{t))), 

where Mk is the operator from Lemma 12.11 associated to the bounded domain K . Then 
div b(t, x) = and b(t, x) = M(t)x + c(t) on dfl for every t > 0. 
If we set u = u — b, then u satisfies (12. ip if and only if u satisfies 

in (0, oo) x Q, 

in (0, oo) x Q, 
on (0, oo) x dn, (2.2) 
for t G (0, oo), 

in Q, 

where f(x) := u (x) — b(0,x) and 

F x (t,x) := Ab(t,x) + (M(t)x+c(t))-Vb(t,x)-M(t)b(t,x)-b(t,x)-Vb(t,x)-b t (t,x). (2.3) 

Note that div / = and that the compatibility assumption ensures that even / G L p a [Q). 
Our approach to system (12. 2p is based on linear operators of the form 

C(t)u(x) = (Aui(x) + (M(t)x + c(t), Vui(x))y ^ - M(t)u(x), t > 0, (2.4) 

and perturbations of the form 

B(t)u(x) = -b{t, x)-Vu-u- V6(t, x), t > 0, 



i*t - Am - (M(t)x + c(t)) ■ Vu + M{t)u 
+b(t, x)-Vu + u- Vb(t, x) + u ■ Vu + Vp 

divM 

u(t, x) 

lim u(t, x) 



Fi{t,x) 







kl 



u(0,x) = f 
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where u = (ui, . . . , Ud) and x is an element from Q, a bounded domain D C M. d or M. d . Note 
that the operators £(•) are of Ornstein-Uhlenbeck type (cf. |12|). In the case of exterior 
domains Q we define the //-realizations of £(•) as 

V(L n (t)) := {n G W /2 ' p (fi) d n W 1,p (fi) d : M(t)x ■ V« G £ p (fi) d }, 
L n (t)u := C(t)u, 
and the perturbed operators are defined by 

V(L a>b (t)) ■= V(La{t)), 
L n , b {t)u := jC(t)u + B(t)u. 

In the following, to simplify our notation, we do not distinguish between L p (Q) d and L P (Q) 
and sometimes write 

L n (t)u(x) := Au(x) + (M(t)x + c(t)) ■ Vu(x) - M(t)u(x), t > 0, x e Q. 

With these linear operators the linearization of the system (12. 2p for some initial time s > 
is now given by 

u t - La jb (t)u + Vp = 0, in (s, oo) x fi, 

divw = 0, in (s, oo) x Q, , , 

it = 0, on (s, oo) x <9Q, ^ ' ' 

u(s, •) = /, in £1 

As usual in the theory of the Navier-Stokes equations we shall later work in the space 
of all solenoidal vector fields in L p . Therefore we set 

V(A n (t)) := {u e W 2 ' p (n) d n Wq' p (QY n LP(fi) : M(t)x ■ V« G 
A n {t)u := P n L n (t)«, 

and 

P(A n , 6 (0) := V(A n (t)), 
A a , b {t)u := ¥ n L njb (t)u. 

By applying the Helmholtz projection to (12. 7p the pressure p can be eliminated and we 
may rewrite the equations as an non-autonomous abstract Cauchy problem 

u'(t) = An, b (t)u(t), 0< s <t, 

u(s) = f. [ ■ } 

It directly follows from j8] that for fixed s > the operator Aq^s) generates a Co- 
semigroup on L P (Q), 1 < p < oo, which is however not analytic. Therefore we cannot 
apply standard generation results for evolution systems of parabolic type (we refer to the 
monographs [TH] and |22] for more information on this matter). Moreover, we note that the 
domain of An^t) depends on time t. Therefore, to overcome this difficulty and in order 
to discuss well-posedness of ( I2.10p we introduce the regularity space 

Y n := {u G W 2 ' p {Q) d n Wl' p (Sl) d n L P {Q) : \x\V Ui (x) G L p {Q) d for i = 1, . . . , d} 

which is contained in D{A^{t)) for every t > 0. 
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Our first main result is the existence of a strongly continuous evolution system on Z£(0), 
1 < p < oo, that solves the Cauchy problem (12.10P on the regularity space Yq. This 
directly implies well-posedness of (12.10p . Moreover, we obtain L p -L q smoothing properties 
and gradient estimates for the evolution system. This is a priori not obvious, since the 
evolution system is not of parabolic type. 

Theorem 2.2. Let O C M. d be an exterior domain with C 1,1 -boundary and 1 < p < oo. 
Then there exists a unique evolution system {T^(t,s)}(t, s )GA on 2v£(0) with the following 
properties. 

(a) For (t, s) G A, the operator Xh,&(i, s) maps Yq into Yq. 

(b) For every f eYq and s > 0, the map t i— > Tn^(t, s)f is differentiable in (s, oo) and 

J^h, 6 (t, s)f = A Q>b (t)T n>b (t, s)f. (2.11) 

(c) For every / G 1q and t > 0, the map s i— > Tn ib (t, s)f is differentiable in [0, £) and 

^T a , b (t, s)f = -T a , b (t, s)A a , b (s)f. (2.12) 

(d) For T > and 1 < p < oo there is a constant C := C(T) > such that for every 

f e 2£(o) 

||VTn, 6 (t,a)/|| p , n < C(t-s)-*||/|| Pi n, (t, s ) g A T . 

(^ej Lei T > and 1 < p < q < oo. TVien i/iere exists a constant C := C(T) > such 
that for every f G 2^(0) 

||2h, 6 (*, s)/||,, n < C(t - S )-*(H) ||/|| p>n , (t, s) G A T , (2.13) 

||VT n , 6 (i, s)/|| g , n < C(t - s )-Ki-iB ll/l^n, s ) g A t . (2.14) 

Remark 2.3. (a) The analogous result holds for the evolution system {T^(£, s)}(t. s )eA 
associated to the operators Aq(-). 
(b) If we denote the evolution system on 2^(0) by {T^ b (£, s)}(£ )S ) 6 a, then the family of 
evolution systems is consistent in the sense that 

T* >b (t, s)f = T^t, s)f, (t, s) G A, / G 2£(0) n 24(0), 

holds for 1 < p, g < oo. 

The basic idea to prove Theorem 12.21 is to study first the whole space case M. d and the 
case of a bounded domain D, which is done in Sections [3] and |4] respectively. Then in 
Section [5] we use some cut-off techniques to construct the evolution system for the exterior 
domain O. Our method, which was already presented in [12] for treating non-autonomous 
Ornstein-Uhlenbeck equations, then also allows to obtain the L p -L q smoothing properties 
and the gradient estimates from the corresponding estimates in M. d and D. 
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Next we come back to the nonlinear problem (12.21) . Again by applying the Helmholtz 
projection Pq we can rewrite this system in abstract form as 

u'(t)-A Q , b (t)u(t)+¥ a (u-Vu)(t) = PnF^t), t>0, 

u{0) = f, 

where Fi is given in ()2.3p . By the Duhamel principle (variation of constant formula) this 
problem can be reduced to the integral equation 

u(t) = T ntb {t,0)f- f T nib (t,s)¥ u (u-Vu)(s)ds (2.16) 
Jo 

T Q , b (t,s)P Q F 1 (s)ds, t>0, 

in L p (R d ). For given < T < oo, we call a function u G C([0, T ); L p (R d )) a mild solution 
of (I2.15P if u satisfies the integral equation (12. 16ft . By adjusting Kato's iteration scheme 
to our situation the existence of a unique local mild solution follows. 

Theorem 2.4. Let 2<d<p<q<oo and f G L P (Q). Then there exists T > and a 
unique mild solution u G C([0,T]; L^(Q)) of \2. 15\) . which has the properties 

t^Xit) E C([0,T};LI(Q)), (2.17) 

ti(?-?)+3 Vu(t) G C([0, T]; L p {Q) dxd ). (2.18) 

If p < q, then in addition 

t%(p~*)\\u(t)\\ q , Q + ^\\Vu(t)\\ p , n ^0 ast^0 + . (2.19) 

3. The linearized problem in R d 

In this section we study the linearized problem in the whole space M. d . This situation 
was already studied in detail by the first author in |11| and by Geissert and the first author 
in [B|. Here we recall the main results. For this purpose we set 

V(L Rd (t)) := {u G W 2 ' p (R d ) d : M(t)x ■ Vu G L p (R d ) d }, 

(3.1) 

L R d(t)u := C{t)u. 

and 

V(A Rd (t)) := V(L Rd (t))nL p (R d ), 
A R d(t)u := C{t)u. 

Note that A R d(t) is indeed an operator on L p (R d ), since an easy computation yields that 
div A R d(t)u = for all u G C^° a (R d ). In particular, the operator L Rd (t) commutes with the 
Helmholtz projection F Rd in R d . Moreover we introduce the regularity space 

Y Rd := {u G W 2 ' p (R d ) d n L p (R d ) : \x\V Ul {x) G L p (R d ) d for i = 1, . . . , d}. 

For every £ > 0, the space Y R d is contained in T)(A Rd (t)). 



THE OSEEN-NAVIER-STOKES FLOW 9 

In the following we denote by {U(t, s)}t, s >o the evolution system in IR d that satisfies 

fZ7(i,s) = -M(t)U(t,s), t>s, 
U(s,s) = Id. 

Now, for / G LP(M. d ) and s > 0, we set T K d(s, s) = Id and for (t, s) G A we define 

T R d(t, s)f(x) = (k(t, s, •) * f)(U{s, t)x + g(t, s)) x G R d , (3.3) 

where 

k(t, s, x) := - 1 f/(t, s )e-i«*"i--->, x G M d , (3.4) 
(47r)2(det g M )2 

g(t,s)= / C/(s, r)c(r)dr and Qt,s = / r)U*(s, r)dr. (3.5) 

For the derivation of this solution formula we refer to [6j Section 3]. The explicit formula 
now allows to prove the following result (see [HUH] and [T2"| Section 2] for details). 

Proposition 3.1. Let 1 < p < oo. Then the family of operators {T R d(t, s)}u,s)eA defined 
in (13. 3p is a strongly continuous evolution system on LP(M, d ) with the following properties. 

(a) For (t, s) G A, the operator T R d(t, s) maps Y R d into Y R d. 

(b) For every f G Y R d and every s G [0, oo), the map t i— > T R d(t, s)f is differentiable in 
(s, oo) and 

J^ R „(t, s)f = A md (t)T md (t, s)f. (3.6) 

(c) For every f G Y R d and t G (0, oo), the map s h-> T Rd (t, s)f is differentiable in [0,t) 
and 

^T Rd (t, s)f = -T Rd (t, s)A Rd (s)f. (3.7) 

(d) Let 1 < p < q < oo. Then there exists a constant C > such that for every 

f e LP(R d ) 

\\T R d(t, s)f\\ q < C(t - s)-i(i-l) ||/|| p , (t, s) G A, (3.8) 

\\VT Rd (t, s)f\\ q < C(t - s )-t(i-§H ||/|| p , ( t , s ) G A. (3.9) 



By Proposition 13.11 and the fact that L m d(t) = A m d(t) on L£(R d ) it follows that the 
solution to the problem 

u t — L R d(t)u = 0, in (s, oo) x M. d , 

divu = 0, in(s,oo)xM d , (3.10) 

u(s,-) = f, inR d , 

for s > 0, is given by u(t,x) = T R d(t, s)f(x). So in the whole space case the pressure is 
constant, thus we may assume that the pressure is even zero. 



(4.1) 



(4.2) 
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4. The linearized problem in bounded domains 
In this section let D C M. d be a bounded domain with C 1,1 -boundary. We set 
V{L Dih {t)) := T>(LD,b) '■= W 2,p (D) d fl Wo ,p (D) d , 
L D)b (t)u := C D)b (t)u. 
Moreover, we define the operators 

V(A D>b (t)) := V(A Djb ):=W 2 *(D) d nWt' p (D) d nI%(D), 
A D , b (t)u := F D L D , b (t)u. 

In a bounded domain also the coefficients of the term M(t)x ■ V are bounded. Thus, it 
follows directly from the classical perturbation theory for the Stokes operator that for fixed 
s > the operator (AD t b{s),T>(Ajj tb )) generates an analytic semigroup (cf. [SJ Proposition 
3.1]). Moreover, by our assumptions on the coefficients we obtain that the map 1 1— > Ao,b(t) 
belongs to C 1 (1R + , ££{J}{An,b)i L p a (D))). Thus, the following result follows from the theory 
of parabolic evolution systems (see (T8| Chapter 6] and |10[ Section 2.3]). 

Proposition 4.1. Let D C M. d be a bounded domain with C 1 ' 1 -boundary and 1 < p < oo. 
Then there is a unique evolution system {T£> b (t, s)}(t jS )eA on L^{D) with the following 
properties. 

(a) For (t,s) G A, the operator T^^it, s) maps L^D) into T>{A£, b ). 

(b) The map t H- Tn jb (t, s) is differentiable in (s, oo) with values in J£(L p a (D)) and 

^T D>b (t, s) = A Dtb (t)T Djb (t, s). (4.3) 

(c) For every f G T>{A£> b ) and t > 0, the map s h- >■ Tu b {t,s)f is differentiable in [0,t) 
and 

^T Dtb (t, s)f = -T Dib (t, s)A Dtb {s)f. (4.4) 

(d) Let T > 0. Then there exists a constant C := C(T) > such that 

\\T D , b (t,s)f\\ p , D <C\\f\\ p>D , (4.5) 

and 

\\T D , b {t,s)fh, P ,D < C(t- S y l \\f\\ PtD . (4.6) 
for all f G LP(D) and all (t, s) G A T . 
For the following first estimate let us recall the Gagliardo-Nierenberg inequality 

\\fh, q ,D < C\\fC^ D \\f\\^ D 9 \ for all / G W"*{D), (4.7) 

where p > 1, q > 1, r > 1, < 9 < 1 and k - | < 9 (m - - (1 - 9)f. 

Corollary 4.2. Let T > and 1 < p < q < oo. Then there exists a constant C := C(T) > 
such that for every f G TP^D) and (t, s) G 

(a) \\T D , b (t,s)f\\ q , D <C(t-sr^)\\f\\ p , D , 
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-11/ 


\\p,D, 


d ( 1 


- 1 ) 1 


2 \p 


q) 2 



p,D- 

Moreover, 

\\Tn, b (t,s)fh, P <C\\f\\ k , p , (t,s)eA T , 
for all f G W k ' p (D) d n I%{D), k = 1, 2, and 

\\T D , b (t, s)f\\ 2>p < C(t - s)-3 (i, s) G At, 

/or a// / G Vr 1 ' p (L>) d n LP(D). 

Proof. Let us assume first that < ^ — ^ < |. Then < | ^~ — M < 1. Hence (a) follows 

by applying (STTJ) with = f (± - ±J , k = 0, m = 2, r = p and (T4~6]) . 

Assume now that 4 < 1 - 1 < > Set - = 1 + Then < 1 - - < \. So, by the first step 

a p q — a r q a p r — a 1 J r 

and ( 14. 7p . we obtain 



\\T D , b (t,s)f\\ q , D = \\T Dib (t,8 + *?)T Dtb (s+¥,s)f\ 



q,D 



_d( 1_ 1 
t - S \ 2 ^ 8 



< cf — Fa* (* + ¥>*)/!,. 



< c 



t — s\ 2 ^ r q > f t — s\ 2 ^ p r > 



2 J V 2 

d ( 1 1 ' 



r,£> 



By iterating this argument we obtain (a) also for ^ — ~ > | . 

Assertion (b) follows from (14. 7p with 6 — ^, p — q — r, k — 1, m — 2 and ( 14. 6p . 
It follows from (a) and (b) that 

\\VT D>b (t,s)f\\ q , D = \\VT D>b (t,s + !f-)T Dib (s + ^,s)f\\ q)D 

_i 

< C (^ir) 2 |l T A6(« + ¥^)/|U 

< c( t _ s )-S(§-*H u/ii^ ( t ,«)eA T) 

and this proves (c). 

For the last assertions we refer, for example, to [T8| Corollary 6.1.8]. □ 



By Proposition 14. II the solution to the problem 



u t — L,D,b(t)u + Vp = 0, in (s, oo) x D, 

divw = 0, in (s, oo) x D, 

m = 0, on (s, oo) x 3D, 

u(s,-) = /, inD, 



(4.8) 
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for s > 0, is given by u(t, x) = To,b(t, s)f(x). From (14. 8 j) and the fact that Ut G L P (D) 
it follows that Vp(t) = (Id — Fo)Lo,b{t)u(t). Since the pressure is only unique up to an 
additive constant, in the case of a bounded domain we can always assume that 

p e L P (D) := ju G L P (D) : J udx = j . 

By the abstract theory of parabolic evolution systems it is clear that u G C 1 ((s, oo); L p (D))r\ 
C((s,oo);T>(Ad^)). So in particular, we can conclude that Vp G C((s, oo), L p (D) d ). By 
Poincare's inequality we can conclude that also p G C((s, oo), L P (D)). 

Lemma 4.3. Let 1 < p < oo and D C ff 1 fc a bounded domain with C 1,1 -boundary. Let 
(u, p) fre the unique solution of (14. 8 j) toii/i p G -^o(-D) and Ze£ 7 G (1 + -, 2). T/ien /or T > s 
there exists a constant C := C(T, 7) > wii/i 

||p(t)IU<C7(t-s)-i||/|U 

/or a// £ G (s, T). 

A similar estimate was proved in [HI Lemma 3.5] for the solution to the corresponding 
resolvent problem (see also [19] for the case of the Stokes operator). However, our proof 
follows the ideas in [21| (see also (2"U| Section 4]). 

Proof of Lemma^^ Let ip G C™(D) and set ip = tp - \D\^ J D pdx. Then, in L p ' (D) 
with ~ + ^7 = 1, there exits a unique solution ip G W 2,p '(-D) of the Neumann problem 

f = <p in A 

[ Vf = on^D, 1 J 

which satisfies the estimate 

for some constant C > (cf. [2"T| Proposition 5.5]). We obtain now 

(p, <p) D = (p, p) D = (p, A^) D = -(Vp, Vip) D = -((Id - W D )L D)h (t)u{t), V^) D 
= -<A«(t), V^)d - ((M(t)z - c(t)) • Vu(i), W>d + (M(t)u(t), V^) D 

+ (b ■ V«(t), W>)r, + ■ V6, V^) D 
= -<Vu(t) ■ 1/, V^) dD + (V«(t), V 2 ^) D - ((M(t)x - c(t)) • V«(t), W>) D 
+ (M{t)u{t), V^) D + (b ■ Vu(t), V^) D + («(*) ■ V6, V^) D . 

By using the embedding Hp +£ ' p (D) m> L p (dD) for < £ < 1 - i (cf [El Theorem 1.5.1.2]) 
we obtain 

|(P,V ? )d| < C(||VM- v\\p,dD + ||Vu|| P) d+ IMkp'.D 
< C||w|| 1+ i +e p D ||y||p',D. 

By taking now e < 1 — - and 7 = 1 + - + e, the assertion follows from Proposition 14. 1[ 
Corollary 14.21 and simple interpolation . □ 
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5. The linearized problem in exterior domains 



In this section we prove Theorem [221 The general idea is to derive the result for exterior 
domains from the corresponding results in the case of R d and bounded domains by some 
cut-off techniques. For this purpose let R > be such that O U supp &(•,£) C B(R) for 
every t > 0. We then set 

D = tin B(R + 8), 

k x = nnB(R + 2), 

K 2 = {x E n : R + 2 < |x| < R + 5}, 

K 3 = {x E n : R + 5 < \x\ < R + 8}. 

We denote by Bj for % E {1,2,3} the operator defined in Lemma 12.11 associated to the 
domain K\. Moreover, by {T R d(t, s)}(t,s)eA we denote the evolution system in L^(M. d ) from 
Proposition 13.11 and by {To, b (t, s)}(t, s )eA the evolution system in L p a (D) for the bounded 
domain D from Proposition 14.11 

Next, we choose cut-off functions 77 E C°°(R ) such that < <p,£,rj < 1 and 

(p(x) := 
t(x) := 



1, 


\x 


1 > i? + 4, 


0, 


\x 


< R + 3, 


1, 1 


X 


>R+1, 


0, 1 


X 


< R, 


1, 


x\ 


<R + 6, 


0, 


x\ 


>R + 7. 



7)(x) := 

I w ) |A| ^ T I ■ 

For a given / E L^(Q) we now define functions f R E L^(R d ) and fo E L%(D), respectively, 
by 

f M ._ / - ii((ve) • /x*), x e n, 

and 

/d(x) := 77(x)/(3:) - B 3 ((V77) • /)(z), x E D. 

With partial integration we see that for every / E V{A^ b {t)) we have f R E V(A R d{t)) and 

f D eV(A D , b (t)). 

Now for (t, s) E A, we define the operator W(t, s) by setting 

W(t, s)f := ipT Rd {t, s)f R + (1 - <p)T D>b (t, s)f D - B 2 ((V<^) • (T Rd (t, s)f R - T D>b {t, s)f D )) 

for every / E L p a (Vt). By Lemma 12.11 it is clear that W(t,s)f E L£(fi). Moreover, for 
/ e V(A Qtb (t)) it follows from the properties of {T Rd (t, s)} (M)eA , {T D>b (t, s)}( M)eA and 

the operator B 2 that W(t,s)/ G H /2 ' p (fi) H W 1,p (n) H Z£(fi) holds. Moreover, a short 
calculation yields 

VW{t, s)f = ipVT Rd (t, s)f R + (1 - <p)VT Djb (t, s)f D 

+ V</? (T Rd (t, s)f R - T D>b (t, s)f D ) 

- V [B 2 ((VyO • (T Rd (t, s)/r - T A6 (t, s)/d))] . 
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Thus, it follows that W(t, s)f G V(A n<b (t)) if f e V(A n , b (t)). 

Let us set Un(t) := T K d(t, s)/r and Uo(t) := To,b{t, s)/d and let prj be the pressure that 
is associated to ud- We may assume that J D p£>dx = 0. 

A short calculation yields that for some initial value / G Yq, the function u[t) := W(t, s)f 
solves the inhomogeneous equation 



Ut - L& >b (t)u + Vp v = -/, in (s, oo) x Q 

divu = 0, in (s 

u = 0, on (, 

it(s,-) = /, infi 



divw = 0, in (s, oo) x f2, (5 1) 

u — 0, on (s, oo) x cftl, 



where Vp^, := V ((1 — </>)p.d) and 

/ := F(t, s)f : = +2 (VT Md (t, s)/ fl - VT D>b (t, s)f D ) ■ (V<p)* 

+ (Ay? + (M(t)s + c(t)) ■ (Vy>)) (T R „(t, s)/r - T A6 (t, s)/ D ) 
+ B 2 ((V^) ■ (ftT R d(t, s)/ fl - dtTi),b(t, s)f D )) (5.2) 
- L n , 6 (t)B 2 ((V^) • (r Rd (t, s)/r - T D b (t, s)f D )) 
+ (V^)pz) 
=: Jx + J 2 + J 3 + J 4 + J 5 . 

Here we use the fact that supp b(t, •) C B(R) for every t > and the expression of 

Certainly, the function F(t, s)/ in (15. 2p is well-defined for every / G an d s ) A. 

Later, we need a certain decay of F(t, s) in t, stated in the following lemma. 

Lemma 5.1. Let 1 < p < oo. For every f G L^(Q) we have 

F(-,-)/GC(A;L^) d ). 

Moreover, let T > fre /ixed and iet 7 G (1 + -, 2). Then 

\\F(t, s)f\\ p , a < C{t - H/IU, (t, s) G A T , (5.3) 
/or some constant C := C(T, 7) > 0. 



Proof. Let us start with the norm estimates for Ji and J 2 . By using Proposition 13.11 and 
Corollary 14.21 we obtain 

\\h\\ P ,n < C(t-8)-t\\f\\ pfl and \\I 2 \\ p>n < C\\f\\ p , a . 

In order to estimate the norm of Z4, let us first note that Lemma |2~T1 implies that Z/Q i b(t)B 2 G 
^(<' p (if 2 ),I/(K 2 )). Thus again by using Proposition 13.11 and Corollary 14.21 we obtain 

ll-^Hp.n < C||/|| Pl n- 
We come now to the term J 3 . At first we note that we can write 
B 2 ((Vy>) • (d t u R (t) - d t u D {t))) 

= B 2 ((V^) • {L Rd {t)u R {t))) - B 2 ((V<^) • (L D , b (t)u D (t) - Vp D )). 
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Now, for a test function ip G C^°(R d ), we have 

\({V<p) ■ L Rd (t)u R (t),^)\ = \(L Rd (t)u R (t),ij(Vcp))\ 

< \(Vu R (t), V U>(V<p)))\ + |<(M(t)a; + c(t)) ■ Vu R (t) - M(t)u R (t), ^(Vcp)}\ 

< C|hii(*)||l,p|| / 0||l 1 p' + C f ||«fl(*)l|lrf»,Jf 2 ||V'llp',Jra 

<C||^(t)|| liP ||V||i, P ', 
where - + -4 — 1. This shows that 

||(V^) •L Rd (t)w il (t)|| w -i, P(Rd) < C\\u R (t)\\ hp 
holds. Analogously, we obtain 

||(V</?) • LD.ftWwDWHw-i.Pp) < C (\\uD(t)\\ ltPtD + ||Pd|| P ,d) • 

Now, Lemma 12.11 together with Proposition 13.11 Corollary 14.21 and Lemma 14.31 yield 

||/3|U<c(t- s )-i||/|u. 

From Lemma [4.31 we can conclude that 

||/5lU<C(t-s)-i||/|U. 

This proves (15. 3p . 

The continuity of the map A 3 (t, s) y F(t, s)f follows from the strong continuity of 
T R d(-, •) and &(■, •), the continuity of the pressure p^(-) together with the properties of 
the operator Bi stated in Lemma 12.11 □ 



Applying to ( 15. ip we have 

Aa,b(t)u - F n f, in (s, oo) x tt, 

(5.4) 





= A Qib (t)u-F n f, 


in 0, 


oo) 


x f2, 


divw 


= o, 


in 0, 


oo) 


x f2, 


u 


= o, 


on (s 


,oo] 


1 x dQ 


u(s, •) 


= /, 


in Q. 







It is clear, that if an evolution system {Tn } b(t, s)}^. s )eA exists on L^(Q), then the solution 
u(t) to the inhomogeneous equation ( 15. 4 p is given by the variation of constant formula 

u(t) = T u ,b(t, s)f - I T Q>b (t, r)F n F{r, s)fdr. (5.5) 

J s 

The integral in (I5.14p exists because of Lemma 15.11 This consideration suggests to consider 
the integral equation 



Ta,b(t, s)f = W(t, s)f + I T a , b (t, r)F n F(r, s)fdr (t, s) e A, f e lf a (0). 

J s 



(5.6) 



Let us now state a lemma which will be very useful in the proof of Theorem 12.21 For the 
proof we refer to |12) . 
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Lemma 5.2. Let X\ and X 2 be two Banach spaces, T > and let R : At —> J£(X 2 , Xi) 
and S : At — > Jzf (X 2 ) be strongly continuous functions. Assume that 

\\R(t, s) U{x 2 , Xl ) < C (t - s) a , \\S(t, s)\\ nx2) < C (t - sf, (t, s) £ At, 

holds for some Co = Cq(T) > and a, (3 > — 1. For f £ X 2 and (t,s) £ At, set 
T (t,s)f := R{t,s)f and 

T n (t, s)f := J T n ^(t, r)S(r, s)/ds, n £ N, (t, s) £ A T . 

TTien i/iere exists a constant C > swc/i t/iat 

00 

X] ||r B (t, < C(t - S ) a ||/|U 2 , (t, «) G A T . (5.7) 

n=0 

Moreover, if a > 0, i/ie convergence of the series in ( f5. 7]) zs uniform on At- 
Proof of Theorem \2.S[ Let T > 0. By using Proposition 13.11 and Corollary 14.21 we have 

s)/|U < C\\f\\ p , a for / £ LP(fi), (t, s) £ A T . 

So, by (15. 3p . we can apply Lemma 15.21 with i? = W, S = FqF, a = 0, = — | and 
Xi = X 2 = L P (Q). Thus, for any / £ L%(Q), the series YlT=o Tk(t, s)f converges uniformly 
in At, where T (t, s)f = W(t, s)f and 

T k+1 (t,s)f = J T k (t,r)F Q F(r : s)fdr, (t,s) £ A T , / £ L p (R d ). (5.8) 
Since T > is arbitrary, 

00 

Tn,6(t,a) := ^T fc (t,a), (t, s) £ A (5.9) 

k=0 

is well-defined. It is easy to check that Tn^t, s) satisfies the integral equation (|5.6j) . 
Moreover, from the strong continuity of W(-, •) and (15. 3p we deduce inductively that 7fc(-, •) 
is strongly continuous and hence, by the uniform convergence of the series we get the strong 
continuity of Th^-, •). 

In order to show that {Tn,&(i, s)}(t, s )eA leaves Yq invariant, we proceed as in the proof 
of Theorem 3.1 in [12] and consider the Banach space 

Z : = {/ £ Wl' p (Sl) d n L P {Q) : |x| ■ Vfi(x) £ L p (fi) d for z = 1, . . . , d} 

endowed with the norm 

\\f\\z:= ||/||i lPl n+|||x|-V/|| Pl n, f^Z. 

Proposition 13. 1\ Corollary 14.21 and (15. 3p permit us to apply Lemma 15.21 with X\ = X 2 = 
Z, R = W, S = F Q F, a = and f3 = -\. So, we obtain that T a>h (t, s)f £ Z for all / £ Z 
and (t, s) £ A. Moreover, by taking X 1 = W 2 > p (Q) d , X 2 = W^ p (n) d n L P (Q), R = W, S = 
PqF, a = — |, (3 = — ^ and applying Proposition 13.11 Corollary 14.21 and (15. 3p . it follows, by 
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Lemma E21 that T a , b (t,s)f G W 2 *{tt) for all / G ^(0)^1^(0) and (t,s) G A. This 
yields that {Th,&(£, s)}(t, s )eA leaves Yq invariant and 



J2l\\Ut,s)fh, P ,n+\\\x\-VT k (t,s)f\\ 



n=0 



< CV(1 + (t - s)-3)(||/|| 1J ,, n + ||M ■ V/|| p ,n), (t, a) G A r , / G Y n . (5.10) 

Let us now prove that for every / G Yq and for every s > fixed, the map £ i— > TQ^(t, s)f 
is differentiable on (s, oo) and that f)2.1ip holds. For / G Yq we compute 

^T (t, s)/ = A n)b (t)T {t, s)f - F n F(t, s)f 



d t T 1 (t,s)f = A sllb (t)T 1 (t,s)f + F a F(t,8)f- f F n F(t,r)F n F(r,s)fd 

J s 

d t T 2 (t, s)f = A Q>b (t)T 2 (t, s)f + J ¥ n F(t, r)F u F(r, s)fdr 



I I F n F{t,r 2 )F n F{r 2 ,r 1 )F n F{r 1 ,s)fdr 2 dr 1 . 
' s J n 



Inductively we see that 

n n 

d t T k (t, s)f = A n>b (t) ]T T k (t, s)f - R n (t, s)f (5.11) 



k=0 k=0 

holds for n G N, where 
-t rt 



Rn(t, s)f:=f I ... I F n F(t, r n )F n F(r n , r B _i) . . . V a F(n, s)f dr n . . . dr 2 d n . 
We estimate now the norm of R n (t, s)f. By Lemma [5.11 we obtain 
\\Ri(t, s)f\\ p ,n < C 2 fit - r)"i(r - S )-idr||/|| p , n = C 2 B(1 - 7 /2, 1 - 7 /2)(t - sf^\\f\\ p ^ 

||ifc(t, s)f\\ p ,a < C 3 B(1 - 7 /2, 1 - 7 /2) f (t - rf^r - *)-idr||/|| ftn 



= C 3 B(1 - 7 /2, 1 - 7 /2)B(l - 7 /2, 2 - j)(t - s) 2 ~'^ lu llp , n . 
Inductively we see that 

\\R n (t, s)f\\ pfl < C n+1 B(1 - 7 /2, 1 - 7 /2)B(l - 7 /2, 2 - 7) . • . 

. . . B(l - 7/2, n - (n 7 )/2)(f - s^^H/IUo 



< \ v ; +1) ; r (*-')- - ( 5 - 12 ) 



C n+1 r(l -j/2) n u (n+l)-, 

[n-^]! 

holds for n G N. Here the constant C may change from line to line. From estimate (I5.12p 
it follows that ||-R„,(t, s)|| p tends to zero as n — > 00. Since we used Lemma [5.21 we know 
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that the convergence of YlT=o Tk{t, s)f is uniform in for any / £ Yn and any T > and 
so, by using (15.10P and the closedness of An sb (t) we can conclude that 

r, OO OO 

- J2 ut, s)f = A n , b (t) m, s )f 

fc=0 k=0 

holds and this proves (12. lip . 

Let us show now the differentiability of the map s t— > Tn,&(t, s)f on [0,t) for t > and 
/ £ Yq. First we note by a short calculation that for / £ X>(Ln &(s)) 

L Rd (s)f R = (Ln, b { 8 )f) R + 2V/ • (VO* + (Af + (M(s)z + c(s)) ■ V$f 

- L Md (s)B 1 ((ve) • /) + Bi((V£) • L n ,6(s)/) 

and 

W*)/d = + 2Vf ■ (V V )* + (At? + (M(s)x + c(s)) • Vt?)/ 

- L A6 ( S )B 3 ((Vr?) • /) + B 3 ((Vt7) ■ LaMf). 

In the following we set Vp := (Id — Pfi)£n,b(s)/. Then we can conclude that for any 
/ £ V(An >b (s)) we have 

F Rd (L n , b (s)f) R = F Rd [CA Q>b (s)f - B^VO^M*)/) + £ V P " Bi((VOVp)] 
= (An, 6 /)ii-PR d [pVe + Bi((VOVp)], 

and analogously 

^D{La,b{s)f) D = (A n , 6 /) D -P^[pV77 + B 3 ((V77)Vp)]. 
Thus, for / £ Yhj we obtain 
d 

— W(t, s)f = -<pT R d(t, s)A Rd (s)f R - (1 - (p)T D)b (t, s)A Djb (s)f D 

+ B 2 ((Vp) • (T Rd (t, s)A R d(s)/.R - T A6 (i, s)A a6 (5)/d)) 
= -^, S )An, 6 (s)/-G(f, S )/ 

where 

G(t, s)f : = <pT R *(t, s)F Rd [2V/ • (VO* + (Af + (M(s)x + c(s)) • VO/ 

-L Md ( S )B 1 ((VO • /) + Bx((VO • W*)/) - pV£ - B^VOVp)] 
+ (1 - y?)T Afe (t, s)¥ D [2V/ ■ (Vr/)* + (At? + (M(s)s + c(s)) ■ Vr/)/ 

- J L A6 ( S )B 3 ((V7 7 ) • /) + B 3 ((Vt7) • Ln fi (s)f) - pVr? - B 3 ((Vr/)Vp)] 
- B 2 {(Vp) • T Rd (t, s)P M , [2V/ • (VO* + (A£ + (M(s)x + c(s)) • VO/ 

-L Rd ( S )B 1 ((VO • /) + Bi((VO ■ L Q , b (s)f) - pV£ - Bi((VOVp)]} 
+ B 2 {(V^) • T Dsb (t, s)¥ D [2V/ • (Vr))* + (A V + (M(s)s + c(s)) ■ Vrj)f 

-L A6 (s)B 3 ((V77) • /) + B 3 ((V77) ■ L n , 6 (a)/) - pVr; - B 3 ((Vr/)Vp)]} 
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We estimate now the norm of p in L P (D) since in the expression of G(t, s)f the supports 
of the functions pV£ and pVr) are subsets of D. Analogously to the proof of Lemma [4.31 
we can show that for any ip G C^°(D) 

\&<p)d\ < CII/HtadIMIpvd 
< C f ||/|| 7i p i n||v?||p' ) D 

and hence 

||p|U < C||/|| 7lPl n 

holds, where 7 G (1 + ~, 2). Now, using the same arguments as in the proof of Lemma [5.11 
we obtain that 

\\G(t,s)f\\ JtP , Q <C(t-syi\\f\\ l!P> n 
for some constant C > and / G H 1,p (Q) d fl L%(Q). Now we apply Lemma 15.21 with 
X 1 = X 2 = H™(Q) d n L p a (n), R = S = G and a = = -\. So, the series 

oo 

V(t,s)f:=Y,Vk(t,s)f, (t, s) G A, 

fc=0 

is well-defined and 

\\V(t, s)/|| 7tP , n < C(t - s)-* ||/|| 7iP ,n, (t, s) G A T , (5.13) 
for / G H 1,p (Q) d fl L P (Q). On the other hand, V(-, •) satisfies the integral equation 

V(t, s)f = G(t, s)f + f V(t, r)G(r, s)fdr, (t, s) e A, / e H™(n). (5.14) 



In particular V(t, •)/ is continuous on [0,£) with respect to the L p -norm for any / G 
H™(Q) d n L£(fi) and t > 0. Now, for / G Z£(fi) and (t, s) G A we set 

f(M)/:=W(t,s)/+ / V(t,r)W(r,s)/dr. 



It follows from the continuity of V(t, -)W(-, s)f on [s,t), Proposition 13. 1\ Proposition 14.11 
and ( I5.13P that the above integral is well-defined for any / G L P (Q). Computing the 
derivative with respect to s G [0, t) yields 

^-f(t,s)f = -W(t,s)L n (s)f -G(t,s)f + V(t,s)f - J V(t,r)W(r,s)L n (s)fdr 

t 

V(t,r)G(r,s)fdr 
= -f(t,s)L n (s)f, 

for any / G Yq, due to (I5.14p . From this equality together with (12.111) and since T^(t, s)Yq C 
Yq, (t, s) G A we can conclude that 

|-(T(*,r)T n , 6 (r, S )/) = 
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holds for all / G Yq. This yields that for / G Yq, the function T(t, r)T^b(r, s)f is constant 
on At and thus, by the density of Yq in L p a (Vt) and by the fact that T > was arbitrary, 
it follows that f(t, s)f = T Qjb {t, s)f holds for all / G (t, s) G A. This proves (12T2|) . 

The uniqueness of the evolution system {Tq.&(£, s)}(£ )S ) £ a follows by a similar argument 
from (12. lip and ( I2.12p (see the proof of Theorem 3.1 in [12] for details). 

The estimate in (d) can be obtained by applying Lemma [5.21 with Xi = W 1,P (Q), X 2 = 
LP(tt), R = W, S = ¥ n F, a = -±,/3 = -f, Proposition O and Corollary H2J 

Finally, the first estimate in (e) follows by applying Lemma T5.2I with X\ = L^(Q), X 2 = 

LP(fi), R = W, S = ¥ Q F, a = -f U - ±V /? = -|, Proposition O and Corollary HJ if 
< -j. The case > § can be obtained by iteration as in the proof of Corollary 

p q — a p q a J 1 J 

14.21 The second estimate follows now from the first one and (d) (see the proof of Corollary 
Q . □ 

To conclude this section let us state the following lemma about the behavior of Th,&(£, s) 
near t = s, which will be needed in the proof of Theorem 12.41 

Lemma 5.3. Let s > and f G L^(JY). Then we have 
(a) for 1 < p < q < oo 

lim (t-sf*(p-«)\\T Q7b (t,s)f\\ q! a = 0, 

t—t-s, t>s 



(b ) for 1 < p < oo 



lim (t-s)3||Vrn l6 (t,s)/|| Pi n = 0. 

t— 5>S, t>S 



Proo/. Let / G < t - s < 1 and choose {f n ) ne ® C C c °° CT (fi) C such that 

limn^oo ||/ — / re || P) n = 0. The triangle inequality together with the L p -L q estimates stated 
in Theorem 12.21 imply that there exist constants Ci,C 2 > such that 

(t - s)i(i-l) ||7h,6(*, a)/|| fl .n 

< (* _ s )f || TfV) (t, s )(/ - / n )||, >n + (t- s)KW) ||T^ b (t, s)/ n || ff ,n 

< dll/ - /„|| Pl n + C 2 (t - s)*(H) ||/ n || ffin , 0<t-«<l,neN. 

Hence, lim(t — s) 1 ||Tq s)/|L n = by letting first £ — > s and then n — > oo. The 
second assertion is proved in a similar way (cf. [TTj Proposition 3.4]). □ 



6. Mild solutions to the nonlinear problem 

In this section we finally come back to the nonlinear problem (12. 2 j) and its abstract 
formulation ( I2.15p . Based on the results proved in Section and an adaptation of the 
Kato iteration procedure (see |17| ) we now prove Theorem 12.41 
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Proof of Theorem\2J\ Let d < p < q, f G L£(fi) and T > 0. For < t < T and fc G N we 
define tti(t) := Tn,&(£, 0)/ and 

u fc+ i(t)=«i(t)- / r nj 6(t,a)Pn(«fc-V« fc )(s)da+ / r n , 6 (t, s^i^s) ds. 
Jo Jo 

Let us prove that, for some T > 0, the sequence (uk)k converges in C([0,T ); L%.(Q)) to a 
mild solution u of fl235|) for T G (0,T ). 

We set /3 := f (j - i) and define 

L fc := L fc (T) := sup ^||u fc (i)|| g ,n, 4 := : = sup ^||Vw fc |U 

te(o,T] *e(o,T] 

and 

M fc := M fe (T) := sup ^||« fc+1 (t)-u fc (t)||, in , M' k := Af£(T) := sup t^\\Vu k+1 {t)-Vu k {t)\\ pM 
te(o,T] te(o,T] 

for fceN. 

It follows from the L p -L q estimates (I2.13P and the boundedness of Pq from L r (Q) d into 
Z£(fi) that 

K+i(*)lkn < IM*)|| 9l n + / \\T a , b (t,s)Fn(u k -Vu k )(s)\\ aQ ds 



Jo 

||r ni 6(t,a)P n F 1 ( a )|| n ds 



< r?L x + C (\t-s)-itt-\)\\u k {s)-Vu k {s)\\ r ^s 

Jo 

+C /ViOOILnd* 
Jo 

< t^Lx + C (*-s)-5(r-?)|| Ufc ( s )|| (I)n ||Vu fc (s)|| p , n ds 

JO 

+C / HFiWILnda 
Jo 

< r fi L x + cr fi L h iJ 1t [ {t-s)-^t p 8 - p -^ds + ct- p f ^11^(8)1^,0 ds, 

Jo Jo 
where - = - + -. Since /3 < | we obtain, by multiplying with t 13 , 

^IK+iWIU < + Ci^Jb + C2T, 

and taking the supremum over t G (0,T] yields 

L k + x <L x + C x L k L' k + C % T, (6.1) 
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where Ci, C2 > are constants independent of k G N but depend on T. Similarly, with 
the gradient L p -L q estimates (12. 14ft we have 

ft 

||Vu fc+ i(*)|| p ,n < t-^L[ + C / (t-s)-5(H)-3|| Ufc ( s ) • Vu fc (s)|| P , n ds 

Jo 

+C [ (t-s)-*||F 1 (s)|| J ,,nds 
Jo 

< Hl; + C / (t-s)~«-?\\u k (s)\\ qj n\\Vu k (s)\\ P! nds (6.2) 
Jo 

+C [ (t-sy^WF^s^nds 
Jo 

and hence, as above, we obtain 

L' k+1 < L[ + C 3 L k L' k + C4T, (6.3) 

where C3, C4 > are constants independent of k G N but depend on T. By setting 
R k := R k {T) := max{L fe , L' k }, it follows from (EE]) and O that i? fc +i < #1 + cii?£ + c 2 T 
for some constants ci, c 2 > 1. By Lemma 15.31 for any e > 0, there is T > such that 
Ri < e for T < To. So, for e < it follows by induction that 

Rk < 4e, fc G N, T < min rf , =: T . 

Thus, the sequences 

{t^t p u k {t)) k and (ti->t3Vu fc (t)) fc (6.4) 

are uniformly bounded in L q a (Vt) and L p (Q) dxd respectively for t < T and all G N. The 
continuity of the maps t 1— > t^U\{t) and t H- taVMi(t) at t = follows from Lemma [5.31 
Hence the continuity of t > t"u k (t) and t t— > Vu k (t) follows from similar calculations. 
We prove now that the sequences in (16.41) are Cauchy sequences. Since 

u k ■ Vu k - u k -i ■ VMfc_i = u k ■ V(u k - u k -i) + (u k - Uk-x) ■ VMfc_i 

and by the same computations as above we obtain 

\\u k+1 (t) - u k (t)\\ q> a < / \\Tn :b (t, s)F a ((u k - Vu k )(s) - (u k -i ■ Vu k -x)(s))\\ n ds 



t 

d 



< C I (t-s) * (||u fc (s)|| ai n||V(u fc (s) - u k -x(s))\\ p ,Q 
+ \\u k (s) -u fc _ 1 (s)|| ffj n||VMfc_ 1 (5)|| p ,n) ds 
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and 



||Vu fc+ i(*)-Vu fc (t)||p,n < f ||Vrn,6(* > s)Pn(K-Vu fc )(s)-(M fc _i-VM fc _i)(s))|| a ds 

Jo 

pt 

< C (t_ s )-5i-5 (|K(s)|| ? ,n||V(u fc (s) - « fc -i(s))|| p ,n 
Jo 



Therefore, 



+ \\u k (s) -Mfc_ 1 (5)|| gi n||Vu fc _i(s)||p j n) ds. 



M k < C 5 (M k _ 1 L k + M k _ 1 L' k _ 1 )<2C 5 R 1 (M , k _ 1 + M k _ 1 ), 
M' k < C e (M k _ 1 L k + M k _ 1 L k _ 1 )<2C 6 R 1 (M k _ 1 + M k _ 1 ), 



and hence, 



(M k + M' k ) < 4e(C , 5 + C , 6 )(M fc _ 1 + M^_ 1 ) 
< \{M k _ x + M' k ^) 

for e < min ^ g^+Cs) ' 8ci) an< ^ < ~^ ' wnere ^5> C*6 > are constants independent of k E 
N but depend on T. Thus, (t \-> t fi u k (t)) k converges to some t ^ t^u{t) E C([0, T], L q a (tt)) 
and [t i — y t^'Vuk(t))k converges to some 1 1— > t^v E C([0, T], L p (Q) dxd ). It is now clear that 
v(t) = Vu(t) and u is a mild solution of (I2.15P on [0, T\. 

The case d < p = q follows by taking (3 = for 5 E (0, 1) in the above computations. 
Here we need ^ + | < 1 to make the integral in (16. 2p convergent. 

Moreover, take any p > d and any r > p and using the same estimates as for ||iifc+i(t)|| 9) n 
above, we obtain 

IK*)IU < \\T(i,b(t,0)f\\ P ,n+ f \\T n , b (t,s)F u (u- Vu)(s)\\ a ds + CT 

Jo 

< C\\f\\ p>Q + C( sup ^||u(t)|| ri n)(su P t*\\Vu(t)\\ p>n ) [ (t - s )-£ a -M ds 
te[o,T] te[o,T] Jo 

+CT 

and hence sup tg [ 0T ] < oo. The continuity can be obtained similarly. Thus, 

uEC([0,T],D>(ti)). 

The case d = p = q can be deduced as in the case d < p = q by taking (5 = for 
5 E (0, 1) and making the some computations but with a q' > d instead of q = d, since if 
q = d then the integral in (16. 2\i could diverge. 

The property (12.191) follows from the construction of the solution u and Lemma 15.31 
Finally the uniqueness follows as in □ 
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